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ABSTRACT 


'nN 

promising  approach  to  failure  modeling,  in  particular  to  developing  failure-time  distri¬ 
butions,  is  discussed.  Under  this  approach,  system  state  or  wear  and  tear  is  modeled  by  an 
appropriately  chosen  random  process,  eg,  a  diffusion  process;  and  the  occurrences  of  fata! 
shocks  are  modeled  by  a  Poisson  process  whose  rate  function  is  state  dependent.  The  system 
is  said  to  fail  when  either  wear  and  tear  accumulates  beyond  an  acceptable  or  safe  level  or 
a  fatal  shock  occurs. 


This  approach  has  significant  merit.  First  it  provides  revealing  new  insights  into  most  of 
the  famous  and  frequently  used  lifetime  distributions  in  reliability  theory.  Moreover,  it 
suggests  intuitively  appealing  ways  for  enhancing  those  standard  models.  Indeed,  this. ap¬ 
proach  provides  a  means  of  representing  the  underlying  dynamics  inherent  in  failure  process¬ 
es.  Reasonable  postulates  for  the  dynamics  of  failure  should  lend  credence  to  prediction  and 
estimation  of  reliability,  maintainability,  and  availability.  In  other  words,  accuracy  of  repre¬ 
sentation  could  lead  to  better,  more  reliable  prediction  of  failure. 
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SECTION  1 
INTRODUCTION 


Modeling  of  failure  times  in  reliability  theory  has  generally  proceeded  by  means  of  ad  hoc, 
albeit  clever,  data  analysis  methods.  In  many  cases,  however,  failure  modeling  has  been  done 
without  any  apparent  attempt  to  capture  the  dynamic  relationship  between  the  system  state 
and  system  reliability.  (Notable  exceptions  include  Gumbel  [10],  Kao  [12],  Birnbaum  and 
Saunders  [5],  and  Bergman  [4];  a  good  discussion  is  provided  in  Mann,  Schafer,  and  Singpur- 
walla  [IS].)  Consequently,  it  has  been  very  hard  to  give  derived  models  a  natural  interpreta¬ 
tion.  Furthermore,  there  is  no  intuitively  appealing  wav  to  fine  tune  a  given  model.  The  lack 
of  an  intuitively  appealing  and  naturally  interpretable  methodology  has  motivated  us  to  study 
reliability  models  that  accurately  reflect  the  dynamic  dependency  of  system  failure  and  decay 
on  the  state  of  the  system. 

Two  relevant  stochastic  models  have  aroused  interest  in  the  theoretical  and  applied  model¬ 
ing  communities;  however,  their  rich  structure  and  potential  for  applications  to  reliability 
have  not  been  fully  exploited.  The  First  of  these  is  the  celebrated  shot-noise  model  (cf.  Rice 
[19]  and  Cox  and  Isham  [7]  for  background  and  further  references).  The  shot-noise  model 
supposes  that  the  system  is  subjected  to  “shots”  or  jolts  according  to  a  Poisson  process.  A  jolt 
may  consist  of  an  internal  component  malfunctioning  or  an  external  “blow”  to  the  system. 
Jolts  induce  stress  on  the  system  when  they  occur.  However,  if  the  system  survives  the  jolt 
it  may  then  recover  to  some  extent.  For  instance,  a  sudden  and  unexpected  surge  of  power 
in  the  circuit  of  a  control  system  may  temporarily  increase  the  likelihood  of  system  failure, 
but  the  overload  itself  decays  rapidly.  For  another  example,  the  mortality  rate  for  persons  who 
have  suffered  a  heart  attack  declines  with  the  elapsed  time  since  the  trauma.  In  this  case,  the 
heart  actually  repairs  itself  to  a  degree.  The  shot-noise  model  is  both  easily  interpretable  and 
analytically  tractable. 

The  other  relevant  model  involves  the  inverse  Gaussian  distribution.  Under  this  approach, 
system  wear  and  tear  is  modeled  by  a  Brownian  motion  that  has  positive  drift.  The  system 
fails  whenever  the  wear  and  tear  reaches  a  certain  critical  threshold.  For  example,  consider 
a  structural  support  subject  to  loadings  that  vary  both  in  terms  of  size  and  in  time  of 
application.  Each  load  causes  microscopic  cracks  to  form  in  the  material.  Eventually  those 
cracks  coalesce  into  a  critical  break  that  causes  the  support  to  fail.  If  the  successive  loadings 
and  their  induced  cracks  are  assumed  to  be  stochastically  independent,  it  is  then  plausible  to 
model  this  process  by  a  Brownian  motion  with  positive  drift.  Under  this  modeling  assumption, 
the  time  to  material  failure  corresponds  to  the  first  passage  time  of  the  candidate  Brownian 
motion  to  the  critical  level,  and  this  first  passage  time  has  an  inverse  Gaussian  distribution 
(cf.  Karlin  and  Taylor  [13,  p.  363]).  Practical  experience  has  shown  that  this  model  provides 
both  good  fit  and  easy  interpretation.  Moreover,  the  distribution  is  extremely  tractable  from 
the  viewpoint  of  statistical  analysis  (cf.  Folks  and  Chhikara  [8]). 


The  rich  structure  of  these  models  suggests  that  an  appropriate  conceptual  framework  for 
reliability  modeling  is  the  following: 

Suppose  that  a  certain  component  in  a  physical  system  begins  operating  with  a  given 
strength  or  a  given  “operational  age"  (eg,  extent  of  wear  and  tear  or  stress)  that  can  be 
measured  in  physical  units.  Suppose  that,  as  time  goes  on,  component  wear  or  stress  builds 
up  (loss  of  strength  with  increasing  age),  perhaps  in  a  random  way.  (The  concept  of  wear  or 
stress  buildup  is  dual  to  that  of  declining  strength.)  For  instance,  consider  the  thickness  of 
tread  on  a  tire  (or  brake  lining)  or  the  level  of  fluid  in  a  hydraulic  system.  The  tread  wears 
down  with  use.  and  there  may  be  gradual  loss  of  fluid  from  the  hydraulic  system.  Assume 
that  this  wear  may  be  offset,  but  only  in  part,  by  maintenance  and  repair.  Such  considerations 
suggest  modeling  component  strength  (or  susceptibility  to  failure)  by  a  stochastic  process  X= 
{X(t),  t  ^  0}  with  starting  state  corresponding  to  the  initial  level  of  strength  (or  initial 
operational  age).  This  process  X  should  tend  to  drift  downward  (decrease)  with  time  as  wear 
builds  up;  if  X  is  the  operational  age  or  wear-and-tear  process  then  it  should  tend  to  drift 
upward.  The  component  may  fail  when  either  wear  alone  has  reduced  strength  below  some 
safe  level  (as  in  structural  material)  or  at  the  epoch  of  occurence  of  some  hazardous  event 
(eg,  an  external  shock)  severe  enough  to  overcome  current  strength.  We  denote  by  r  the  time 
of  passage  of  the  X  process  to  the  critical  level.  For  the  examples  cited  above,  hazardous 
events  might  be  the  tire’s  abrupt  encounter  with  a  sharp  portion  of  road  surface  or  the  rupture 
of  tubing  in  the  hydraulic  system  due  to  an  external  blow.  (It  is  clear  from  these  examples 
that  the  rate  of  fatal  shocks  should  be  modeled  as  a  decreasing  function  of  component  strength 
or  an  increasing  function  of  component  wear  or  stress.)  We  denote  by  k(x)  the  Poisson  killing 
rate  associated  with  state  x,  and  by  T  the  time  to  failure  of  the  component.  With  the  above 
conventions  and  modeling  assumptions  in  place,  we  can  express  the  probability  of  surviving 
beyond  time  t,  starting  with  strength  or  operational  age  x,  as  follows: 

(1)  PX(T>t)  =  EX(expj-^  k(X(s))ds|  I|r>t|  ] 

The  next  step  in  the  modeling  process  is  the  selection  of  a  mathematical  model  for  the 
strength  (wear)  process.  Two  viable  candidate  classes  are  diffusion  processes  and  shot-noise 
processes.  We  first  explore  the  diffusion  model  in  some  detail. 

Suppose  that  component  strength  evolves  in  accordance  with  a  diffusion  process  X  =  {X(t), 
t  >  0}  having  drift  parameter  n(x)  and  diffusion  coefficient  <r2(x)  in  state  x  >  0.  (A  compre¬ 
hensive  introduction  to  the  subject  of  diffusion  processes  is  provided  in  Chapter  15  of  Karlin 
and  Taylor  (14).)  Then-u(x)  can  be  interpreted  as  the  rate  at  which  wear  builds  up  in  state 
x.  Alternatively,  if  component  wear  builds  up  according  to  the  diffusion  process  X  then-M(x) 
can  be  interpreted  as  the  rate  at  which  strength  declines  in  state  x.  If  T  is  the  time  to  failure 
of  the  component,  we  assume  that  P(T  <  h  I  X(s)  =  x)  *  k(x)h  +  o(h)  for  each  time  point 
s,  so  that  k(x)  can  be  interpreted  as  the  Poisson  rate  of  occurrence  of  a  traumatic  shock  of 
magnitude  sufficient  to  overcome  or  “kill”  a  component  of  strength  x. 


Now  let  w(x.t)  be  the  probability  that  a  component  of  strength  x  survives  beyond  time  t. 
ie.  wtx.t)  =  P*  (T  >  t).  Indeed,  w(x.t)  coincides  with  the  right  side  of  equation  ( l)  if  r  is  the 
first  passage  time  of  the  diffusion  X  to  the  critical  level.  It  follows  from  the  backward 
differential  equation  for  the  Kac  functional  of  the  diffusion  process  X  (cf.  [14.  pp.  222-224]) 
that  w(x.t)  satisfies 

3  w  i  x.t)  . .  dwix.t)  j:(x)  d*w  (x.t) 

,2)  - =  -k( x)  w  (x.t)  +  u(x)  -  +  -  — ^ - 

3t  3x  2  3  x- 

The  initial  condition  for  this  differential  equation  is  determined  by  the  critical  strength  or 
wear  threshold  A.  We  allow  0  <  A  <  4-  ®.  If  the  diffusion  represents  strength  then  w(x,0) 
=  1  if  x  >  A  and  0  otherwise.  If  the  diffusion  represents  wear  or  stress  then  w(x,0)  =*  l  if  x 
<  A  and  0  otherwise.  Examples  of  subclasses  that  lead  to  mathematically  tractable  solutions 
of  the  backward  differential  equation  (2)  for  the  state-dependent  survivor  function  w(x,t)  are 
presented  in  Table  l-l .  In  the  table  the  diffusion  process  X  represents  wear  or  stress  buildup. 
These  subclasses  overlap  one  another.  Moreover,  we  certainly  do  not  imply  that  (2)  has  a 
tractable  solution  for  every  choice  of  infinitesimal  parameters,  killing  rate  function,  and 
failure  threshold  in  these  subclasses.  From  a  failure  modeling  standpoint  there  is,  rather,  a 
rich  and  Interesting  variety  of  cases  for  which  explicit  solutions  can  be  Identified. 
Furthermore,  the  widespread  availability  of  excellent  and  efficient  computational 
algorithms  to  solve  (2)  approximately,  makes  it  feasible  to  experiment  with 
various  choices  of  Infinitesimal  parameters  for  modeling  purposes. 


Table  1-1.  Model  Subclasses  and  Parameters 


Subclass 

InfinitMimal 

Parameters  uix),  (Pix) 

Killing 

Function  k(x) 

Deterministic 

Wear 

a2  (x)aO 

arbitrary 

Constant  Killing 
Rate 

arbitrary 

k  (x)aX 

Infinite  Level 
of  Wear  and  Tear 

arbitrary 

arbitrary 

Failure 
Threshold  A 


arbitrary 


SECTION  2 

THE  DETERMINISTIC  SYSTEM  STATE  PROCESS 


For  this  modei  the  buildup  of  wear  (or  loss  of  strength)  is  assumed  to  be  deterministic,  and 
thus  <r2(x)  ■  0.  Under  this  assumption  the  system  state  process  X  =  |X(t),  t  2=  0}  satisfies 
the  following  equation: 

(3)  X(t)  3  (  m<X<s)Ms  +  * 

*•0 

Alternatively,  X  satisfies  the  differential  equation 

(4)  dX(t)  =  M(X(t))dt 

with  initial  condition  X(0)  =  x.  The  probability  of  system  survival  through  time  t  is  then  given 
by 

(5)  P(T>t)  *  exp  I-  fWfsJidsl  I 

Equation  (5)  can  be  interpreted  in  terms  of  standard  reliability  theory.  To  simplify  the 
notation  suppose  for  now  that  A  =*  +  ®  and  x  =  0,  and  write  P(T  >  t)  in  place  of  w(x,t).  First 
we  need  to  recall  the  failure  rate  function  r(t),  which  is  usually  defined  by 

(6)  r(t)  =-(dP(T  >  t)/dt)/P(T  >  t) 

Integrating  r  over  [0,t]  and  assuming  that  P(T*0)  *  0  shows  that 

^  r('s)ds  =  -Ln(P(T>t)) 
o 

Consequently 

(8)  exp|-  y  r(s)ds|=  P(T>t) 

o 

Comparing  (5)  and  (8)  we  see  that 

(9)  r(t)  =  k(X(t)) 

Therefore  the  failure  rate  at  time  t  is  equal  to  the  killing  rate  that  corresponds  to  the  “system 
state"  at  time  t. 


Moreover,  equations  (4)  and  (9)  enable  us  to  reinterpret  standard  reliability  distributions 
in  terms  of  a  killing  rate  function  and  a  system  evolution  process.  To  illustrate  this,  we  will 
now  give  new  interpretations  to  some  well  known  distributions  in  reliability  modeling. 

Makeham’s  Life  Distribution 

(10)  P(T  >  t)  =  exp  {-bt -a(exp(ct)- 1 )} 

Then 

(11)  k(X(t))  =  r(t)  =  d(bt  4*  a(exp(ct)-l))/dt  =  b  +  acexp(ct). 

By  taking  the  killing  rate  proportional  to  the  stress,  eg. 

(12)  k(x)  =  x 
we  find  that 

(13)  X(t)  =  b  +  acexp(ct) 

Equation  (4)  implies  that 

( 14)  n(b  +  acexp(ct))  =  ac2  exp(ct) 

Thus 

(15)  M(x)  -c(x-b) 

Interpret  the  “stress”  function  m(x)  as  saying  that  an  individual  is  born  b  operational  units 
old  and  declines  at  a  rate  proportional  to  chronological  age.  If  we  set  b=0  in  (10)  we  obtain 
the  Gompertz  life  distribution.  Gompertz  [9]  described  the  conceptual  basis  for  his  distribu¬ 
tion  as  assuming  “...the  average  exhaustion  of  a  man’s  power  to  avoid  death  to  be  such  that 
at  the  end  of  equal  infinitely  small  intervals  of  time  he  lost  equal  portions  of  his  power  to 
oppose  destruction  which  he  had  at  the  commencement  of  these  intervals.”  Gompertz  thus 
assumes  that  an  individual’s  operational  age  is  proportional  to  his  chronological  age.  Make- 
ham’s  improvement  to  the  Gompertz  distribution  takes  into  account  the  fact  that  a  person 
is  bom  b  operational  units  old. 

Weibull  Distribution 

(16)  P(T  >  t)  =  exp(-atb) 

Then 

(17)  k(X(t))  =  r(t)  =  d(atb)/dt  =  abtb"' 

Suppose  that  b  1.  Taking  k(x)  =  xb'\  gives 

(18)  X(t)  =  t(ab)1,<b"1>. 

Setting  c  =  (ab),/(b'0  will  simplify  (18)  to 

(19)  X(t)  -ct 

This  implies  the  wear  rate  function  is  constant,  ie. 


(20)  4(x)  ■  c 

Ifb  =  l  then  k(X(t))  =  a  tor  all  t.  Thus,  either  the  killing  rate  function  is  constant  or  the 
wear  rate  is  zero  or  both.  Adopting  the  exponential  distribution  as  a  failure  model  is  a  strong 
assumption  indeed! 

Rayleigh  Distribution 

(.21)  P(T  >  t)  =  exp{-(at  +  bt2/ 2)} 

Then 

(22)  k(X(t))  =  r(t)  =  a  +  bt 
Taking 

(23)  k(x)  =  x 
gives 

(24)  X(t)  =  a  4-  bt 
and  thus 

(25)  m(x)  =  b  and  x<,  =  a 

Gumbel  Distribution 

(26)  P(T  >  t)  »  1  -  exp(-exp(-(t-a)/b)) 

Then 

(27)  exp  [  -  ^  k  (X(s))dsj  =  1  -  exp  (-exp(-(t-a)/b)) 

o 

or 

(  28)  ^  k(X(s))ds  =  -  Ln  (  1  -  exp(-expHt-a)/b))> 
o 

Hence 

(29)  k(X(t))  =  {exp(-exp(-(t-a)/b))exp(-(t-a)/b)/b|/{!-exp(-exp(-(t-a)/b))} 

*  Ein(exp(-(t-a)/b))/b 

where  Ein(x)  is  the  Einstein-Planck  function  (cf.  Abramowitz  and  Stegun  [  1  j )  defined  by 

(30)  Ein(x)  *  x/(exp(x)  -  l) 
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Setting 

(31)  k(x)  =  Ein(x)/b 
we  ;ee  that 

(32)  X(t)  =  exp(-(t-a)/b) 
and  consequently 

(33)  M(X(t))  =»-exp(-(t-a)/b)/b 

where  n(x)  —  -x/b.  In  contrast  to  the  Makeham,  Weibull,  and  Rayleigh  examples,  it  is 
appropriate  to  view  this  X  process  for  the  Gumbel  distribution  as  representing  loss  of  strength. 

Repeating  the  foregoing  development  with  A  <  +  «  (ie.  a  finite  critical  failure  threshold) 
yields  truncated  versions  of  these  well-known  distributions. 

Thus,  within  the  framework  of  our  general  approach,  a  “first  order  model"  (taking  ^(x) 
■  0)  not  only  yields  some  classic  failure  distributions  but  provides  new  and  revealing 
insights  into  their  structure.  It  Is  Important  to  note  that  while  the 

decompositions  between  state  process  and  killing  rate  function  presented  above 
are  somewhat  arbitrary,  this  should  not  be  the  case  with  an  actual  application.  In 
many  reliability  modeling  efforts  there  will  be  a  natural  candidate  for  system 
state,  and  the  corresponding  state  process  will  be  observable.  This  observability 
permits  the  study  of  conditional  probability  of  failure  given  the  state  process. 

Indeed,  the  conditional  distribution  of  failure  due  to  trauma  Is  simply  the 
distribution  or  the  first  event  In  a  nonhomogeneous  Poisson  process  whose  rate  at 
time  t  Is  equal  to  k(X(t)).  It  Is  this  fact  that  allows  the  analyst  to  decompose  the 
parameter  estimation  problem  Into  two  distinct  and  largely  Independent 
statistical  estimation  problems;  namely,  the  estimation  of  the  killing  rate  function 
and  the  estimation  of  the  system  state  process  parameters. 

Furthermore,  this  decomposition  has  significant  statistical  implications.  In 
most  Instances  the  state  of  the  system  at  death  will  contain  all  the  relevant 
Information  about  the  killing  rate  function.  Independent  of  the  time  of  death. 
Consequently,  post-mortem  analysis  will  yield  most  If  not  all  of  the  relevant 
Information.  The  actual  time  of  death  becomes  a  less  significant  factor. 

We  will  consider  these  Important  Implications  In  a  future  paper. 
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SECTION  3 

NONDETERMINISTIC  WEAR  (STRENGTH)  MODELS 
(WHY  DO  WE  NEED  THEM?) 


The  deterministic  wear  model  presented  above  is  satisfactory  when  used  in  either  predict¬ 
ing  the  behavior  of  very  large  ensembles  of  units  (such  as  assessing  the  price  of  warranties 
for  mass  produced  items)  or  in  predicting  the  behavior  of  units  whose  wear  process  is  nearly 
constant.  In  many  cases,  however,  neither  of  these  conditions  is  satisfied.  For  example,  the 
setting  of  safety  standards  for  the  replacement  of  airplane  tires  is  a  case  in  point.  The  wear 
of  tires  is  obviously  quite  random;  tire  wear  depends  on  many  factors  such  as  weather 
conditions,  runway  surface,  airplane  load,  etc.  The  determination  of  the  optimal  tire  replace¬ 
ment  policy  requires  a  reliability  model  for  the  tires.  Inaccurate  reliability  analysis  will  result 
in  either  unnecessary  tire  replacement  or  in  tire  failures  (with  possibly  catastrophic  conse¬ 
quences).  The  utility  of  a  ’‘dynamic"  reliability  model  is  that  it  allows  the  safety  engineer  to 
define  a  replacement  policy  in  terms  of  tire  condition  rather  than  in  terms  of  operational  age. 
By  exploiting  the  information  contained  in  tire  condition,  it  is  possible  to  simultaneously 
minimize  the  average  number  of  replacements  while  decreasing  the  risk  of  failure. 

One  plausible  way  to  obtain  reliability  models  that  satisfactorily  capture  the  dynamic 
connection  of  component  reliability  with  component  condition  is  to  model  wear  or  stress  by 
a  diffusion  process.  A  diffusion  process  can  be  conceived  as  a  continuous  approximation  of 
a  first  order  stochastic  difference  equation,  ie,  something  of  the  form 

(34)  X(n+I)  -  X(n)  =  M(X(n))  +  <x(X(n))Z„ 

where  {Z„,  n  ^  1 }  is  a  sequence  of  independently  and  identically  distributed  random  variables. 

Interpret  (34)  as  saying  that  during  one  unit  of  time  (and  starting  from  level  x)  the  wear 
will  increase  on  average  n(x)  and  that  the  standard  deviation  from  this  average  increase  in 
wear  is  <r(x).  If  <r(x)  is  identically  zero,  then  (34)  reduces  to  a  first  order  difference  equation 
whose  continuous  approximation  is  given  by 

(35)  X’(t)  =  M(X(t)) 

If  we  assume  that  the  sequence  (Zn,  n  >  1}  is  normally  distributed,  then  (34)  is  simply  a 
discrete  approximation  to  a  diffusion  process.  It  is  certainly  plausible  to  model  the  distribution 
of  the  {Zn,  n  >  1}  sequence  by  some  other  distribution  (eg,  a  gamma);  however,  the  normality 
assumption  seems  to  produce  mathematically  more  tractable  equations.  Furthermore,  there 
is  an  enormous  literature  on  diffusion  processes  which  can  be  utilized. 

Several  subclasses  of  our  diffusion  approach  in  failure  modeling  seem  to  yield  tractable 
problems.  They  are  discussed  in  the  following  paragraphs. 

First  Passages  With  Constant  Killing 

Consider  a  unit  that  is  replaced  when  its  strength  dips  below  a  certain  level  or  when  it  fails 
due  to  a  catastrophic  shock.  Assume  that  the  probability  of  death  due  to  a  catastrophe  is 


independent  of  the  unit's  condition.  (For  example,  a  car  headlight  is  replaced  as  scon  as  its 
luminescence  drops  below  a  certain  critical  level  or  it  is  broken  by  some  road  shock.)  In  this 
case,  the  model  has  constant  killing  rate  function,  ie.  k(x)  »  X.  Thus  the  probability  that  the 
“component  is  alive  and  well”  at  time  t  is  given  by 

(36 )  PX(T  >  t)  =  EX [I  |  r>t^  exp  (-X  t  )1 

where  r  is  the  time  of  first  passage  of  the  diffusion  process  to  the  critical  level  A.  This 
expression  simplifies  to 

(37)  F‘X(T  >  t)  =  exp(-Xt)  PX(r  >  t) 

Consequently,  it  suffices  to  derive  the  distribution  of  the  first  passage  time  r.  Calculating  the 
first  passage  time  distribution  for  an  arbitrary  diffusion  is  generally  quite  difficult.  In  fact, 
it  corresponds  to  solving  a  second  order  partial  differential  equation  with  nonconstant  coeffi- 
cientsj  in  oarticular,  equation  (2)  with  k(x)  »  0  and  the  appropriate  initial  condition. 

However,  there  Is  a  significant  class  of  diffusions  that  give  rise  to  mathematically 
tractable  first-passage  problems. 


(to  page  3-3) 
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Generalized  Inverse  Gaussian 

Any  generalized  inverse  Gaussian  distribution  with  a  nonpositive  power  parameter  is  a  first 
passage  distribution  for  a  certain  diffusion  process  with  drift  (cf.  Barndorff-Nielsen  et  al  [3]). 
The  generalized  inverse  Gaussian  distribution  has  density 

,  a/ 2 

(c/b)  a_t  -1  i 

(38)  fit)  =  - t  exp<-(bt  +  ct)/2  /  it  >01 

2Ka(v/"Sc) 

Here  K.a(  )  stands  for  the  modified  Bessel  function  of  the  third  kind  (cf.  [1])  with  index  a, 
and  a  is  called  the  power  parameter.  The  domain  of  variation  of  the  parameters  (a.b.c)  is 

f30)  a  >  0,  b  ^  0.  c  >  0 

a  =  0,  b>0.  c>0 

a  =5  0,  b>0,  c  £  0 

This  nch  class  of  distributions  includes  the  inverse  Gaussian  (a  =-1/2),  the  hyperbola 
distribution  (a  =  0),  and  as  limit  cases,  the  gamma  (a  >  0  and  b  =  0)  and  the  reciprocal 
gamma  (a  <  0  and  c  =  0).  (For  a  recent  survey  of  the  main  results  and  related  references, 
see  Jorgensen  [11].) 

For  a  first  passage  example  (cf.  [3])  and  a  striking  Illustration  of  the  richness  of 
this  class,  consider  the  following.  Let  g  be  the  an  arbitrary  strictly  Increasing 
function  on  [O,00)  with  continuous  second  derivative  for  which  g(0)  =  0.  Let 
X(t)  =  g(crB(t)— jxt)  where  a  >  0  and  p.  >  0.  Now  set 

(40)  <t(x)  =  o-g'(g_1(x)), 

and 

(42)  |i(x)  =  <7(x)[~t  +  ^2(Sl] 
ct  2dx 
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Consider  the  diffusion  on  [0,°°)  with  Infinitesimal  parameters  p.(x)  and  (T(x)  so 
defined.  Let  x0  >  0  be  the  initial  position  of  the  process.  Then  the  first  hitting 
time  of  level  0  has  the  distribution  obtained  from  (38)  with  parameter  set 
(a,b,c)  =  (-1/2,  02(xo),  v?/ <J 2). 

For  Instance,  take  g(x)  =  x2.  Then  o^x)  —  4ct2x,  02(x)  =  x/cr2,  and 
jjl(x)  =  —2\ixU2  +  a2'  For  another  example,  take  g(x)  =*  exp(x)  —  1;  this  Is 
geometric  Brownian  Motion  centered  at  the  origin.  Then 

o r(x)  =  o^x  +  l)2,  02(x)  =  f  ancj  y,(x)  =  (— p,  +  cr/2)(x  +  l). 

<T 

(The  transformation  of  diffusion  processes  Is  discussed  In  Karlln[l4].) 


Failure  Due  to  Shock  Only 

In  this  class  of  models,  systems  can  age  indefinitely  (A  =*  ®);  that  is,  a  system  is  not  subject 
to  retirement  but  must  die  in  the  line  of  duty.  As  an  example,  consider  home  appliances,  which 
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are  literally  worked  to  death  at  ages  well  beyond  any  nominal  design  lifetime.  Another 
example  is  automobile  usage  in  third  world  countries,  where  cars  are  driven  500.000  miles 
or  more  before  retirement. 

Mathematically  speaking,  the  probability  of  survival  for  this  class  corresponds  to  the  Kac 
functional  of  the  state  process.  More  specifically,  let  X  be  the  state  process  and  let  Y(t)  be 
defined  as  follows: 


-s: 


kl  X(s))ds 


Then  the  probability  of  survival  beyond  time  t  is  given  by 
(-H)  w(x,t)  =  E*[exp(-Y(t))] 

For  a  remarkable  example,  consider  Brownian  motion  with  quadratic  killing. 

To  motivate  this  model  In  a  reliability  context.  Imagine  that  a  system  ages 
(wears)  on  the  average  p.  age  units  per  unit  of  time,  but  that  the  Increase  in  age 
per  unit  time  deviates  from  Its  expected  value  by  a  age  units.  Furthermore, 
imagine  that  the  system  Is  subject  to  externally  generated  shocks,  some  of  which 
are  potentially  fatal.  For  example,  consider  the  human  cardiovascular  system  and 
suppose  that  we  are  measuring  age  In  terms  of  some  operational  criterion  (eg, 
systolic  blood  pressure).  The  human  heart  Is  constantly  subject  to  shocks 
(emotional  distress,  sudden  physical  exertions,  etc).  The  likelihood  that  a  given 
shock  will  cause  heart  failure  depends  on  the  condition  of  the  heart  (and  the 
overall  health  of  the  person),  so  It  Is  reasonable  to  postulate  that  the  occurrence 
rate  of  fatal  blows  Increases  with  age.  Moreover,  the  occurrence  rate  of  fatal 
blows  grows  at  first  very  slowly  with  age,  then  more  rapidly  with  advancing  age. 
The  above  considerations  would  suggest  using  a  convex  killing  rate  function  that 
has  a  minimum  at  0.  A  plausible  candidate  function  Is  k(x)  —  Xx2  . 


Let  X  be  Brownian  motion  with  drift  p<  and  positive  variance  o-*  ,  and  Initial 
position  x  ,  le, 


t-10)  X(t)  ~  <rB(t)  +  jxt  +  x 
Let  the  killing  function  be  given  by 

(46)  k(x)  =  Ax2 

The  functional  Y(t)  is  then 

2 

(47)  Y(t)  =  X  V  (X(S)1  ds 

o 

For  this  important  case,  the  probability  of  survival  up  to  time  t  (starting  in  state  x)  is  given 
by 

(48)  w(x,t)  =rsech  (at v^rl  exp  /  -  +  _JfL_  tanh(aty3C)  - 

L  J  '  2a2 


+  sech  (oty/TX)  -  tanh(at>/5X)  | 


where 


sech  y  *  2/(exp(y)  +  exp(-y)] 
and 

tanh  y  =  [exp(y)  -  exp(-y)]/[exp(v)  +  exp(-y)] 

The  calculation  of  this  expression  for  w(x,t)  is  a  nontrivial  exercise.  To  obtain  the  formula 
given  by  the  right-side  of  (48)  we  first  expressed  X(t)  by  means  of  the  Karhunen-Loeve 
expansion  (cf.  Ash  and  Gardner  [2]).  The  calculation  then  required  the  use  of  special  function 
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theory  and  the  calculus  of  residues  as  it  applies  to  the  summation  of  series.  Finally,  we  checked 
our  calculation  by  verifying  that  w(x.t)  satisfies  the  backward  differential  equation  (2)  with 
k(x)  =  Xx2.  m(x)  *  <r(x)  *  <r,  and  initial  condition  w(x.O)  =  l.  The  derivation  of  (48)  and 

some  properties  of  the  distribution  are  presented  in  [15]. 

Shot  Noise  Model  for  System  Stress 

We  continue  in  the  vein  of  the  preceding  example  (failure  due  to  externally  induced  trauma 
only)  but  return  to  the  formalism  of  equation  (1).  Suppose  that  system  wear  or  operational 
age  is  modeled  by  a  shot  noise  process  X  =  (X(t),  t  ^  0}.  For  example,  suppose  the  system 
is  subjected  to  shots  or  jolts  according  to  Poisson  process  with  rate  X.  Suppose  that  if  a  jolt 
of  magnitude  D  occurs  at  epoch  S  then  at  time  S  4- 1  the  contribution  of  the  jolt  to  the  system 
stress  is  Dh(t),  where  h  is  a  nonnegative  function,  vanishing  on  (-®,0),  which  tends  to  0 
sufficiently  fast  as  t  -*  ®  for  several  integrals  of  powers  of  h.  over  (0,«).  to  be  finite.  In  other 
words,  shot-induced  stress  is  additive  and  decays  with  time  according  to  the  rate  function  h 
(recall  the  cardiac  example  mentioned  earlier).  Thus,  if  (Sn,  n  >  1}  are  the  epochs  of  shot 
occurrences  and  f  D„,  n  ^  1 }  the  magnitudes  of  the  successive  jolts,  then  the  “residual”  system 
stress  at  time,  say  X(t),  is  given  by 

OO 

(49)  X(t)  D„h(t-V 

n=l 


An  intuitively  appealing  and  customary  choice  for  h  is  h(y)  =  exp(-ay)  if  y  >  0,  and  h(y) 
=  0  if  y  <  0,  ie,  exponential  decay.  Suppose  that  system  failure  is  proportional  to  system 
stress,  say  k(X(t))  =  X(t)  in  equation  (1),  and  X(0)  =  0.  Then  the  time  to  system  failure  T 
is  the  epoch  of  the  first  count  in  a  doubly  stochastic  Poisson  process  with  rate  function  ?X(t), 
t  >  0}.  In  particular,  if  the  shot  epochs  and  magnitudes  are  independent  and  the  shot 
magnitudes  are  mutually  independent  with  a  common  distribution  having  Laplace  transform 
<t> ,  then 

(50)  P(T>t)  =  E  [exp  |  X(s)ds|  ] 


=  exp(  -  X  t)  exp 


jx  J*  ^Jyh(z)dz^dyj 


This  distribution  has  failure  rate 
(31)  X  [l-0(  fh(y)dy)] 

The  formulation  is  intuitively  appealing  while  leading  to  tractable  results.  Indeed,  system 
stress  X(t)  has  a  limit  distribution  as  t  -*>  ®.  In  particular,  if  D  has  an  exponential  distribution 
with  parameter  n  and  h(y)  =  exp(-ay)  for  y  >  0,  then  the  limit  distribution  of  X(t)  is  gamma 
with  location  parameter  \i  and  shape  parameter  X/a.  This  is  a  remarkable  result  in  that  a 
gamma  distribution  with  arbitrary  shape  parameter  appears  as  a  limit  distribution  in  a 
physical  model  (For  background  and  further  discussion  on  the  structure  of  shot- 
noise  distributions,  see  BondessonfeU 
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Various  enhancements  are  possible,  including  alternative  choices  for  the  attenuation  or 
recovery  function  h  and  modeling  the  pattern  of  shot  occurrences  by  a  nonhomogeneous 
Poisson  process,  a  renewal  process,  a  semi-Markov  process,  or  a  cluster  point  process.  Implica¬ 
tions  of  the  shot-noise  formulations  will  be  explored  In  a  future  naner. 
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